rs={'all',''};%文件前缀
for rno=1:size(rs,2)
    sname=rs{rno};
datay=load([sname,'datay.txt']);
datax=load([sname,'datax.txt']);

theta=datay(:,2)*180/pi;
phi=datay(:,3)*180/pi;
rindex = boundary(theta,phi); 
figure
pps=polyshape(theta(rindex),phi(rindex));
plot(pps)
hold on
plot(theta(rindex),phi(rindex))
fig=figure
scatter(theta,phi)
hold on
plot(theta(rindex),phi(rindex))
if strcmp(sname,'')
    saveas(fig,'realboundnocons.svg')
end
end